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Abstract 

The question of the "physical meaning" and "origin" of the Bose-Einstein (BE) factor 
in the fiuctuation-dissipation theorem (FDT) is often raised and this term is sometimes 
interpreted as originating from a real harmonic oscillator composition of the physical 
system. Such an interpretation, however, is not really founded. Inspired by the famous 
work of Caldeira and Leggett, we have been able to show that, whenever linear response 
theory is applicable, which is the main hypothesis under which the FDT is established, 
any generic bosonic and/or fermionic system at temperature T can be mapped onto a 
fictitious system of harmonic oscillators so that the suscettivity and the mean square of 
the fiuctuating observable of the real system coincide with the corresponding quantities 
of the fictitious one. We claim that it is in this sense, and only in this sense, that the BE 
factor can be interpreted in terms of harmonic oscillators, no other physical meaning can 
be superimposed to it. At the best of our knowledge, this is the first time that such a 
mapping is explicitly worked out. 



The fluctuation-dissipation theorem (FDT) [I] is very general and applies to a broad 
variety of different physical phenomena. An often raised question concerns the physical 
meaning and/or origin of the Bose-Einstein (BE) distribution factor which appears in the 
relation between the power spectrum of the fluctuating quantity and the corresponding 
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generalized suscettivity. Sometimes this term is interpreted as due to an harmonic oscil- 
lator composition of the considered physical system. Such an interpretation, however, is 
not supported by the derivation of the theorem itself. Moreover, the FDT applies to any 
generic bosonic and/or fermionic system. 

Far from being an academic question, this issue is of very practical importance in 
many different contexts ranging from solid state physics to astrophysics and cosmology. 
Actually, from a real understanding of the origin of this term often depends the correct 
physical interpretation of theoretical and experimental results (see below). In this letter 
we solve this important interpretation issue. 

To this end, inspired by seminal papers of Caldeira and Leggett [21 [3], in the following 
we shall be able to establish a new very general result (a mapping between any generic 
system at temperature T and a fictitious system of harmonic oscillators), whose interest 
goes beyond the specific application to the FDT investigated in the present work. 

Given a generic bosonic or fermionic system which interacts with an external field f{t) 
through the interaction term V = —f{t) A, where A is an observable (a bosonic operator) 
of the system, the fluctuation-dissipation theorem states that, whenever linear response 
theory is applicable, the mean square of the Fourier transform A{uj) of A(t) is related to the 
imaginary part x'^ {^) of the corresponding (Fourier transformed) generalized suscettivity 
by the relation [1] : 

(AH:.)) = %:(-)coth (^^) =2nxM (1 + , (1) 

where (3 = l/kT, with T the temperature of the system and k the Boltzmann constant. 

Let us consider, for instance, a resistively shunted Josephson junction [4j and, more 
specifically, the power spectrum Si{uj) of the noise current (fluctuation) in the resistive 
shunt (dissipation). When applied to this case, the theorem takes the form {R is the 
shunt resistance) [5]: 

The power spectrum Sj{uj) has been measured [4] and good agreement between the ex- 
perimental results and Eq. (I2|) has been found. The specialized literature [H [5l [6l [7] often 
presents the ^ term in parenthesis as due to zero point energies and the experimental re- 
sults |lj as a measurement of zero point energies. In fact, the term in parenthesis in Eq. ((21) 
coincides with the mean energy of an harmonic oscillator of frequency u in a thermal bath. 
The same holds true for the general case of Eq. ([1]), where the term in parenthesis is the 
mean energy of an harmonic oscillator in Hw units, i.e. the BE distribution function. 

The question that we would like to answer is if it is physically founded to interpret this 
BE term as coming from an harmonic oscillator structure of the system. In particular, 
we would like to understand if the agreement between the experimental results [4] and 
Eq. ([2]) can be considered as a signature of measurement of zero point energies. 

We anticipate now the results of our analysis. On the one hand, we shall be able 
to establish a mapping between the physical system and a flctitious system of harmonic 
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oscillators in such a manner that the mean square of A{uj) and the related (imaginary 
part of the) generalized suscettivity X^{^) of the (real) system are precisely reproduced 
by considering the corresponding quantities of the fictitious one. At the same time, our 
analysis will clarify that it is only in this sense that the BE factor can be interpreted in 
terms of harmonic oscillators. Therefore, no physical oscillator degrees of freedom of the 
system are involved in Eq. ([I]) and no zero point energies have been measured in [4]. 

Let us begin by briefly reviewing the derivation of the FDT. Consider a macroscopic 
system with unperturbed Hamiltonian Hq under the influence of the perturbation 

V = -f{t)A{t), (3) 

where A{t) is an observable (a bosonic operator) of the system and f(t) an external gener- 
alized forc^. Let \En) be the Hq eigenstates (with eigenvalues En) and {En\A{t)\En) = 0. 
Within the framework of linear response theory, the quantum-statistical average {A(t)) j 
of the observable A(t) in the presence of V is given by 

{Mt))f= f dt'xAt-t')f{t') (4) 

where Xa ~ ^0 the generalized suscettivity, 

xAt - t') = 'fit - t'){[A{t),A{t')]) = -^Gn{t - t') , (5) 

with (...) = E„f?n(K|...|K), Qn = e-^^'-fZ , Z = Ene"^""", GR{t-t') being the 
retarded Green's function and A{t) = Q^Hot/hj^^-iHot/h^ 
Deflning the correlators (from now on t =0): 

G>(t) = (i(t) i(0)) and G<(t) = (i(0) A{t)) , (6) 

so that Gi?(t) = —i6{t){G>{t) — G^{t))^ and the corresponding Fourier transforms, G^{uj) 
and G^{uj) respectively, it is a matter of few lines to show that: 

g>(^) = - ^_^„^^ ImGfi(cu) ; G<(u;) = e-^^G>(u;). (7) 
Finally, by noting that 

(^'M) = ^(G'>(a.) + G<(a;)) (8) 
and that the Fourier transform of Xa(^) Xa(^) — X^l"^) + ^X^'(^) = ~\Gr{'^) we get: 

(•4^(-)> = = fixiMcoth (^1^) = 2fix:(-) (i + ^^) (9) 

^More generally, we could consider a local observable and a local generalized force, in which case 
we would have V = — j fA{r) f (f, t) , and successively define a local suscettivity (see 
Eq. ^ below). As this would add nothing to our argument, we shall restrict ourselves to r-independent 
quantities. The extension to include local operators is immediate. 
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which is Eq. (P), the celebrated FDT. 

As observed by Kubo et al. [8j (and shown in the derivation sketched above), the BE 
factor is simply due to a peculiar combination of Boltzmann factors in Eq. ([9|) and there 
is no reference to physical harmonic oscillators of the system whatsoever. Despite such an 

authoritative remark, some people insist in interpreting the + ^i3hi_i ^ term as related 

to an harmonic oscillator structure of the physical systenjfl. 

To begin with our analysis, it is useful to show that from Eqs. (E]) and (JH]) we can easily 
derive the following expressions for Xa{^) ^ind {A^^u)) : 



h 



(10) 



and 
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h 



Ej — Ej 

h 



(11) 



where Aij = {Ei\A\Ej) . 

In fact, by inserting in Eq. ((5|) the expressions: 

e{t -t') = - /_+^ t^-^ , / = m{E.\ and A{t) = e^Hot/n^^-^Hot/n 



oo 2ni ui+irj ' 

1 r duj e-^-^^*-*') 



h J 2tx uj + if] 



{E,-E,){t-t')/h _ ^-i{Ei-E,){t-t')/h 



(12) 



Then, by defining the Bohr frequencies Uij = {Ei — Ej)/h and making use of the relation 



lim„_>o -^-^ 



P (^) — mb{uj) we have for xj(t — t'): 



x'At-t') = I / — y ^?,|A,,f fe-*(----)(*-*')5(c..)-e-^(-+'^-)(*-^^ 

h J 271 ^ — ^ V / 

= Ij E ^^^i^^^f (^(^ + ^^^0 - ^(^ + ^^•^)) (13) 

which immediately gives Eq. (fTOl) . 

As for Eq. ([TT]) . by inserting the identity in Eq. ([6|) for G>(t) we find: 

G>(t) = {Ait)A{0)) = J^p,{E,\ei^'Ae-i^'\E^){E^\A\E,) = 

id 

= 5^Ae-*(^^-^»)*|A,/. (14) 

id 

^In the case of the measured [4] power spectrum of Eq. ([2]), some authors [51 [TO] interpret this term as 
due to the electromagnetic field in the resistive shunt and therefore the first term in parenthesis of Eq. ^ 
as originating from zero point energies of this electromagnetic field. 
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Working out the similar expression for G'<(t), for the correlation function G{t) we finally 
have: 



G{t) = ^(G>(t) + G<(t)) = i5^p,|A,,p(e-*(^-^^)* + e-*(^-^^)*) 



and taking the Fourier transform: 

G{lj) = 7r^pi\Aij\^ 



Ej — Ei \ r 

^ ' +uj]+6 



As 



G{t) 



h 



+°° rl, I 



Ej — Ej 

h 



and (A^) = G{0), G{lj) is nothing but the spectral density of (A^) 



(15) 



(16) 



(17) 



(18) 



From Eqs. (fT6i ). (flTl) and (fTHl ). we immediately get Eq. (fTTl ). 
For our purposes, it is useful to write Eq. (flTl ) in a different manner. It is not difficult to 
show that we can actually write Eq. ( JTTll aq 



J6. 



{uJij + Lu) + 6 {ujji + CU)] 



TT 



coth ^^^^ ^lil^iiP [S {uJij + u) - 6 {uji + u)] . 



(19) 



After few additional steps we finally have: 

{A\u;)) = TT 



^{Qi - gj)\Aij\'^ coth ( 
j>i ^ 



TT 



[5 {iu - LUji) + 6{iu + ujji)] (20) 
coth ^{Qi - [S {uJ - Uji) -5{uj + Uji)] . (21) 



Similarly, it is easy to see that Eq. ([T0|) can be written as: 



(22) 



j>i 



Starting from Eqs. (l20l) and (l22l) and taking inspiration from previous work of Caldeira 
and Leggett |3j, we shall be able to establish a formal mapping between the real system 



^Obviously, comparing Eq. ifTQ]) with Eq. ifTO)) . we find, as we should, the FDT theorem, but this is not 
our goal. 



considered so far and a system of fictitious harmonic oscillators. A similar mapping, 
restricted however to the T = temperature case, was considered in |3j, where it was also 
noted that the T 7^ case needs separate discussion. The mapping that we are going 
construct in the present work actually deals with the T 7^ general case. At the best of 
our knowledge, this is the first time that such a mapping is explicitly worked out. 

To prepare the basis for the construction of this mapping, let us consider first a real 
system Sosc of harmonic oscillators (each of which is labeled below by the double index 
{ji} for reasons that will become clear in the following) whose free Hamiltonian is: 

j>i ^ J* / 

where ujji are the proper frequencies of the individual harmonic oscillators and Mji their 
masses. Let \nji) {riji = 0, 1, 2, ...) be the occupation number states of the {ji} oscillator 
out of which the Fock space of Sosc is built up. Let us consider also Sosc in interaction 
with an external system through the one-particle operator: 

Vosc = -fit)Aosc, (24) 

with 

Asc = ^ (ttji Qji) ■ (25) 

j>i 

Obviously, the FDT applied to Sosc gives: 

(^L(^)) = %1(^) coth (^^^ , (26) 

but this is not what matters to us. 

What is important for our purposes is that, differently from any other generic system, 
for Sosc we can exactly compute {AI^^{lj)) and Xosd^) from Eqs. ffTOll and ffTTl) because 
we can explicitly compute the matrix elements of Aosc- 

In fact, if we apply Eqs. (fT4l ). ( fTSi l and (fT6l ) to Sosc and replace the double index 
notation introduced above {uji ; riji ; Mji ; etc.) with the more convenient (for the time 
being) and self explanatory one index notation (tui, ^2, ... ; rii, n2, ... ; Mi, M2, ... ; 
etc.), for {Alg^{uj)) we have: 

■)\{ni,n2, ..\A 

osc |mi,m2, ..;| 

ni,n2,.- mi,m2,-- 

x[5{uj + hcui + /2C^2 H ) + 5{uj - liui - I2UJ2 )] (27) 

where 4 = nfe - m^, g^, = e-/3("fc+V2)?!^fe/^^^ = e-^("'=+^/2);i^fc (j^^^g ^^^^ ^-^at in 
this one index notation Aosc is written as Aosc = J2k {(^k(jk) )• Now, as 



{nk\qk\mk) = 
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we immediately get: 



ni,n2,.. mi,m2,.. 



X 



X [5 (u; + liui + /2t^2 H ) + 5{uj - huji - I2UJ2 )] • (29) 



Let us concentrate our attention to the square in the second line of Eq. ( l29l) . Due 
to the presence of the Kronecker deltas, all the crossed terms in this square, i.e. all the 
terms with different values of the index /c, vanish. In other words, the square of the sum 
is equal to the sum of the squares: 

. k ' ^ h^k 



For the same reason, the same holds true for each value of the index k, i.e.: 

2 



(30) 



h 



, 2 



2MkUJk 

h 



hj^k J 



h^k 



Therefore, as h = — n^'k, for (A^^^(u;)) we get: 



(^L(^)) = ^ 5^ (f?„,^?„,---) 5Z«fcl^^f;^(2nfc + l)((5(u;-L^fc) + 5(u; + u;fc)) (32) 

ni,n2,.- fc 



Finally, as Qn^ = 1, the above expression becomes: 

{AlJuj)) = TTY,(^l^^^i^i^-^k) + 6iu + Uk))Y(^nA'2nk + l) (33) 

k rik 

= ^ E (5(c. - .;,) + 5{u + c,)) . (34) 

k 

Going back to the original double index notation, the above equation is written as: 

{AIM) = nY,c.l^J^,^oth(^y6{u;-u,,) + 6{u; + u;,.)) (35) 

= TT coth (^^^ Yl "i* 2M-U- ^^^^ " ^^'^ -S{uj + uji)) . (36) 
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We have just seen that given a real system Sosc of harmonic oscillators and the one 
particle operator Aosc of Eq. (l25|l . for such an operator is possible to evaluate explicitly 
{Al^^{uj)). We find that each of the individual harmonic oscillators gives rise to a term 

coth ( j which in turn comes from the term ^„ . Qnji (2 riji + 1) of Eq. (l33l) . 



X. 



Let us now consider Xosci^)i which (see Eqs. (fTOl) and (129]) ) is nothing but: 



[UJ 



ni,n2,- 



mi ,m2 ,. 



ttfc 



h 



2MkUJk 



X [(5 (CJ + /iCJi + I2UJ2 + ■■■)— 5 {uJ — liUi — I2UJ2 



(37) 



Apart from the factor 1/h, Eq. (l37l) differs from Eq. fl29l) because it contains the difference 
(rather than the sum) of delta functions in the last line. 

If we proceed for Xosd^) have just done for {Al^^{uj)), we immediately note that 

the only difference with the previous computation is due to this minus sign. In fact, its 
presence causes that rather than the combination (2 rifc + 1) of Eq. (I32i l. which comes from 
the sum {rik + 1) + of Eq. (l3T]) . we get the combination (n^ + 1) — = 1. Therefore, 
for x'osci'^) not get the sum (2 Uk + l) = coth (^^y^) of Eq. (l33|) . but rather 



1. Then: 



^ \^ 2 



2MjiUJji 



{S{UJ - LOji) - 6{UJ + Uji)) . 



(38) 



Naturally, comparing Eq. (l36ll with Eq. (l38! ) we see that for Sosc the FDT holds true, as 
it should. However, what is important for our purposes is to note that for this system we 
have been able to compute separately {Alg^{uj)) and x'osd^) found that the coth (^^) 



of each of the 



factor of the FDT originates from the individual contributions coth 

harmonic oscillators of Sosc- 

We are now in the position to build up our mapping. Let us consider the original 
system S, described by the unperturbed Hamiltonian Hq , in interaction with an external 
field f{t) through the interaction term V = — fit) A (see Eq. ([3])), and construct a fictitious 
system of harmonic oscillators Sosc, described by the free Hamiltonian Hose of Eq. ([23l) . in 
interaction with the same external field f(t) through the interaction term Vosc of Eq. ( l24ll , 
with Aosc given by Eq. ( I25l) . where for aji we choose: 



a 



2MjiUjji 

h 



and for the proper frequencies a;,j of the oscillators: 



UJ 



{Ej - E,)/h > 



(39) 



(40) 



8 



with Ei the eigenvalues of the Hamiltonian Hq of the real system. 

Comparing Eq. (l38l) with Eq. (l22|) and Eq. ([36l) with Eq. (j2Tl) . it is immediate to see 
that with the above choices of aji and ujji we have: 

(A'iuj)) = {AlJu;)) (41) 
x'lM = XosM- (42) 

Eqs. (j39l) and ([40l) are the central results of our analysis. These are the equations that 
allow to establish the desired mapping. Actually, with such a choice of the a's and the 
tu's, we have been able to derive Eqs. ([4T]) and ([42l) . What we have just found is that it is 
possible to map the real system S onto a fictitious system of harmonic oscillators Sosc , 

S ^ Sosc, (43) 

in such a manner that x'^i^j) and (A^(u;)) of the real system, i.e. the imaginary part of 
the (Fourier transformed) generalized suscettivity and the mean square of the (Fourier 
transformed) interaction operator A(t), are equivalently obtained by computing the cor- 
responding quantities of the fictitious one. 

It is worth to point out that the key ingredient to construct such a mapping is the 
hypothesis that linear response theory is applicable, which is the main hypothesis under 
which the FDT is established. When linear response theory is not valid, Eq.Q cannot 
be derived. As a consequence, we do not arrive to Eqs. (ITOi ) and (ITTi ) which are crucial to 
build up our mapping. 

Moreover, by considering the "equivalent" harmonic oscillators system Sosc rather than 
the real one, we are somehow allowed to regard the BE distribution factor coth ('^) of 

the FDT in Eq. ([T]) as originating from the individual contributions coth j of each 

of the oscillators of the fictitious system (see above, Eqs. fl35ll . (l36l) and fl38l) ). In this 
sense, such a mapping allows for an oscillator interpretation of the BE term in the FDT. 

At the same time, however, our result shows that this BE factor does not describe 
the physics of the system, i.e. it does not encode any real, physical, harmonic oscillator 
structure of the system (see also the considerations below). In this respect, we have to 
stress that what we have implemented is not a canonical transformation, i.e. it is not a 
transformation which allows to describe the system in terms of new degrees of freedom 
(such as normal modes), but just a formal mapping, a mathematical construct, which can 
be established, we repeat ourselves, only within the framework of linear response theory. 

In our opinion, our finding provides a definite answer to the often raised questions 
of the "physical meaning" or "physical origin" of the BE term in the FDT or, stated 
differently, to the question of whether this BE distribution factor possibly describes the 
physical nature of the system or not [TT] . 

In fact, from the derivation of the FDT, we know that the BE factor derives from 
a peculiar combination of Boltzmann factors (see [8] and Eq. ([9]) above). At the same 
time, we have shown that, regardless the bosonic or fermionic nature of the (real) system 
S, it is always possible to establish a mapping which relates 5 to a system of harmonic 
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oscillators Sosc so that this BE factor can be regarded as "originating" from the individual 
oscillators of the "equivalent" system Sosc- Therefore, it is not the physical nature of the 
system which is encoded in this BE term but rather a fundamental quantum property of 
any bosonic and/or fermionic system: whenever linear response theory is applicable, any 
generic system is, at least with respect to the FDT, equivalent (in the sense defined above) 
to a system of quantum harmonic oscillators. 

Let us summarize now our results. We have found that, when linear response theory 
is valid, for any generic bosonic and/or fermionic physical system S whose interaction 
with an external field f{t) is given by a term of the kind V = —f{t)A, where A is an 
observable of the system, it is always possible to find an "equivalent" system of harmonic 
oscillators such that Xj^i^) and (A^(a;)) of the real system can be obtained by comput- 
ing the corresponding quantities of the fictitious one. The operator that represents the 
physical observable A is such a mapping is the one-particle operator Aosc of Eq. (l25|) with 
the aji and the Uji given by Eqs. (l39l) and (l40ll . In passing (i.e. in order to establish such 
a result), we have shown that for any system of harmonic oscillators interacting with an 
external field through an interaction term of the kind given in Eqs. (I24l) and ( l25l) . the 
BE coth (^) factor originates from the contribution of each of the individual harmonic 
oscillators of the system. As the physical operator A in the mapping is represented by the 
operator Aosc of Eq. (I25ll . we concluded that the BE distribution factor which appears in 
the FDT Eq. ^ can be regraded as originating from the individual harmonic oscillators 
of the equivalent system. Our results clearly indicate that it is only in this sense that this 
BE factor can be interpreted in terms of harmonic oscillators and that no other physical 
meaning can be superimposed to it. 

We believe that our mapping has a broader range of applicability than the worked 
case of the FDT discussed in this letter. Work is in progress in this direction. 
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